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Abstract: The energy evolution of average multiplicities and multiplicity fluctuations in jets produced in 
heavy-ion collisions is investigated from a toy QCD-inspired model. In this model, we use modified split- 
ting functions accounting for medium-enhanced radiation of gluons by a fast parton which propagates 
through the quark gluon plasma. The leading contribution of the standard production of soft hadrons 
is enhanced by a factor \/TV^ while next-to-leading order (NLO) corrections are suppressed by 1/ 
where the parameter A'^s > 1 accounts for the induced-soft gluons in the medium. Our results for such 
global observables are cross-checked and compared with their limits in the vacuum. 
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Recent experiments at the Relativistic Heavy Ion Collider (RHIC) have established a phenomenon of 
strong high-transverse momentum hadron suppression [1], which supports the picture that hard partons 
going through dense matter suffer a significant energy loss prior to hadronization in the vacuum (for 
recent review see [2]). 

Predictions concerning multi-particle production in nucleus-nucleus collisions can be carried out by us- 
ing a toy QCD-inspired model introduced by Borghini and Wiedemann in [3]; it allows for analyti- 
cal computations and may capture some important features of a more complete QCD description. In 
this model, the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) splitting functions q ^ gq and 
9 ^ 99 [4] of the QCD evolution equations were distorted so that the role of soft emissions was en- 
hanced by multiplying the infra-red singular terms by the medium factor Ng. The model [3] was further 
discussed and used on the description of final states hadrons produced in heavy-ion collisions [5]. 

Within the model, we make predictions for the medium-modified average multiplicity Na in quark and 
gluon jets {A = q,g) produced in such reactions, for the ratio r = Ng/Nq and finally for the second 
multiplicity correlators {Na{Na — 1)) /N\, which determines the width of the multiplicity distribution. 

The starting point of our analysis is the NLO or Modified-Leading-Logarithmic- Approximation (MLLA) 
master evolution equation for the generating functional [4] which determine the jet properties at all 
energies together with the initial conditions at threshold at small x, where x is the fraction of the outgoing 
jet energy carried away by a single gluon. Their solutions with medium-modified splitting functions can 
be resummed in powers of y/ag/Ng and the leading contribution can be represented as an exponential 
of the medium-modified anomalous dimension which takes into account the A^s-dependence: 



Within this logic, the leading double logarithmic approximation (DLA, 0{y/T^J^a^)), which resums both 
soft and colUnear gluons, and NLO (MLLA, 0{as)), which resums hard coUinear partons and accounts 
for the running of the coupling constant a^, are complete. The choice dY = d@/@, where B ^ 1 is 
the angle between outgoing couples of partons in independent partonic emissions, follows from Angular 
Ordering (AO) in intra-jet cascades [4]. In order to obtain the hadronic spectra, we advocate for the 
Local Parton Hadron Duality (LPHD) hypothesis [6]: global and differential partonic observables can be 
normalized to the corresponding hadronic observables via a certain constant K. that can be fitted to the 
data, i.e. N^,^ = K. x Ng^q. 

The evolution of a jet of energy E and half-opening angle involves the DLA anomalous dimension 70 
related to the coupling constant through 7q = 2Ncas/iT, with = 2t: /ANc[5q{Y + A), where Y = 
ln((5/(5o) {Q = is the hardness or maximum transverse momentum of the jet), A = ln((5o/A^^^) 
is a parameter associated with hadronization (Qo is the collinear cut-off parameter, kx > Qo, and A 

QCD 

is the intrinsic QCD scale) and /?o = 4^ (t-^c ~ |-^^)' where Tr = nj/2, rij being the number 
of active flavors. At MLLA, as a consequence of angular ordering in parton cascading, the average 
multiplicity inside a gluon and a quark jet, Ng^g, obey the system of two-coupled evolution equations [7] 
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where 7med(«s) can be expressed as a power series of ^/a^/N^ in the symbolic form: 
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which follow from the MLLA master evolution equation for the generating functional; Ng^g = Ng^qiY), 
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which accounts for parton energy loss in the medium by enhancing the singular terms like <^ ^ Ng/x 
X <C 1 as proposed in the Borghini-Wiedemann model [3]. Thus, when Ng increases the DLA becomes 
dominant and energy-momentum conservation plays a less important role. 

For y » Inx ~ ln(l — x), Ng^q{x) (Ng^q{l — x)) can be replaced by Ng^g in the hard partonic splitting 
region x ^ 1 — x 1 (non-singular or regular parts of the splitting functions), while the dependence 
at small x ^ 1 is kept in the singular term <I>(x) Ns/x as done in the vacuum. Furthermore, the 



integration over x can be replaced by the integration over Y{x) 
approximate system of two-coupled equations, 

d^ 



In 



2Ee\ 

Qo J 



dY^^^ 
d^ 



Cp 2 



. Thus, one is left with the 
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with the initial conditions at threshold A^a(O) = 1 and A^^(O) = and the hard constants: 
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ai = 3/4. 



The quantum corrections oc ai , oi in (I5I6I ) arise from the integration over the regular part of the splitting 
functions, they are 0{yJas/Ns) suppressed and partially account for energy conservation as happens in 
the vacuum. 

These equations can be solved by applying the inverse Mellin transform to the self-contained gluonic 
equation which leads to 



N^'{Y) ~ /C X 



ujpo exp 



uj{Y + \) + 
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where the contour C lies to the right of all singularities of Ng{u}) in the complex plane. Since we are 
concerned with the asymptotic solution of the equation as 1" » 1 {EQ » Qo), that is the high-energy 
Umit, the inverse Melhn transform (|7]l can be estimated by the steepest descent method. Indeed, the large 
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Figure 1: MLLA ^ medium-modified average multiplicity as a function of Q = EQ in tlie vacuum 
(Ng = 1) and in the medium {Ng = 1.6 and Ng = 1.8) fornj = 3. 



parameter is Y and tfie function in the exponent presents a saddle point at uq = \/ Ng/ I3q{Y + A), such 
that the asymptotic solution reads 
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N^{Y)c^Kx{Y + \) ^oexpJ^(Y + X), 
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where ui = ^ — The constant o"i is A^^ -independent because itresums vacuum corrections. There- 
fore, the production of soft gluons in the medium becomes exp 



2{Vn;-1)V{Y + X)/Po\ higher 
than the standard production of soft gluons in the vacuum [4]. From ([B and ([8]l one obtains the medium- 
modified MLLA anomalous dimension 7n,ed = '^^w' ~ V'^To ~ '''iTo' which is nothing but the 
MLLA rate of multi-particle production with respect to the evolution-time variable Y in the dense 
medium. In Fig.[T] we display the medium-modified average multiplicity ([8]) with predictions in the 
vacuum {Ng = 1) in the range 10 < Q(GeV) < 200; we set Qq = A^^^ = 0.23 GeV in the hm- 
iting spectrum approximation [7], and take /C = 0.2 from [7]. The values Ng = 1.6 and Ng = 1.8 
in the medium may be realistic for RHIC and LHC phenomenology [3, 5]; the jet energy subrange 
10 < (5(GeV) < 50 displayed in Fig. [U has been recently considered by the STAR collaboration, which 
reported the first measurements of charged hadrons and particle-identified fragmentation functions from 
p-i-p collisions [8] at ^snn = 200 GeV. Finally, the whole jet energy range in the same figure, in particu- 
lar for those values at Q > 50 GeV, will be reached at the LHC, i.e Q = 100 GeV is an accessible value 
in this experiment (see [3] and references therein). 

We find, as expected, that the production of soft hadrons increases as A'^s > 1: the available phase 
space for the production of harder coUinear hadrons is restricted as the model itself states. The medium- 
modified MLLA gluon to quark average multiplicity ratio, r = Ng/Nq = Ng /N^, following from ^ 
and (O reads 
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where we introduced the coefficient ri = ai — ai in the term suppressed by 'yo/y/N^ as Ng > 1. 
Therefore, if compared with its behavior at = 1, we check, as expected from the model [3], that r 
becomes closer to its asymptotic DLA limit ro = Nc/Cp = 9/4, as depicted in Fig.|2l Setting A^^ = 1 
in one recovers the appropriate limits in the vacuum [4, 9, 10]. Finally, the gluon jets are still more 
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Figure 2: MLLA ratio r (O as a function of Q = E@ in the vacuum {Ns 
{Ns = 1.6 and Ns = 1.8) for n/ = 3. 



1) and in the medium 



active than the quark jets in producing secondary particles and the shape of the curves are roughly the 
same. 

The normalized second multiplicity correlator A2 = {Na{Na — 1))/Nj^ defines the width of the multi- 
plicity distribution and is related to its dispersion by the formula D\ = {A2 — 1)N^ + Na [9]. These 
moments, which are less inclusive than the average multiplicity, prove to be /C-independent and therefore 
provide a pure test of multiparticle production. The medium-modified system of two-coupled evolution 
equations for this observable follows from the MLLA master equation for the azimuthally averaged gen- 
erating functional [4] and can be written in the convenient form 
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which proves to be more suitable for obtaining analytical solutions in the following. We use a new 
method to compute solutions at MLLA by replacing N^^ = A2N'^ on both sides of the expanded 
equations at x ~ 1 — x ~ 1. The notations in (IIOII II) follow the same logic as those in (I5I6I) . Applying 
the analysis that led to the system (I5I6I) . we obtain from (IIOIIII) 
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Figure 3: MLLA second multiplicity correlator inside a gluon jet (fTSl ) as a function of Q = EQ in the 
vacuum (Ng = 1) and in the medium (Ng = 1.6 and A'^, = 1.8 for = 3. 



The constant A^,; only affects the leading double logarithmic term of the equations. The terms propor- 
tional to ai, (ai — bi) and ai are hard vacuum corrections, which partially account for energy conserva- 
tion, indeed ~ V^lo the relative correction to DLA is 0{^J'ar[N~s). 
Setting n'^^ = G2Ng in ([T2l ) and making use of dSjl, the system can be solved iteratively by taking terms 
up to 0{as) into consideration. The analytical solution reads. 



while its expansion in the form 1 + 70/ y/N^ leads to 



where the linear combination of color factors can be written in the form 



c, = J-a-4^ + H^^^-^^^). (16) 



We use dEl) and © and substitute A^f ^ = Q2N^ into ([T3] ) such that the solution reads 

where we obtain the combination of color factors 

1 f55 , ATrCf STnCj. 



4nX9^9N,N, 3N,N^)- ^^^^ 

Setting = 1 in ([T5] ) and ([TT] ) we get a perfect agreement with the vacuum results [9]. In Fig.[3]and 
Fig. 151 we compare our results for the medium-modified second multiplicity correlators ([T5] ) and ([TT] ) 
with predictions in the vacuum (Ng = 1) [9] in the limiting spectrum approximation inside the typical 
range 10 < (5(GeV) < 200 for RHIC and LHC phenomenology. Similarly to the MLLA ratio r{Ns), 
Eq. (O, the hard corrections 0(70) are suppressed by a factor l/\/Ns. As expected from the model. 
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Figure 4: MLLA second multiplicity correlator inside a quark jet ([TT] ). 



we check that these results approach their DLA limits when Ng increases; moreover, the multiplicity 
fluctuations of individual events must be larger for quark jets as compared to gluon jets just like in the 
vacuum [9]. Another interesting feature of these observables concerns the shape of the curves. They 
are roughly identical and prove not to depend on the medium parameter Ns. Moreover, there exists 
evidence for a flattening of the slopes as the hardness of the jet Q = EQ increases for Ns > 1 (vacuum 
and medium). This kind of scaling behavior is known as the Koba-Nielsen-Olsen (KNO) scaling [11]: 
it was discovered by Polyakov in quantum field theory [12] and experimentally confirmed by e^e~ 
measurements [13] for the second and higher order multiplicity correlators. 

In this paper we have dealt with the medium-modified average multiplicity and the medium-modified 
second multiplicity correlator in quark and gluon jets at RHIC and LHC energy scales. The starting point 
of our calculations is based on the Borghini -Wiedemann work [3], which models parton energy loss in 
a nuclear medium. The average multiplicity is found to be enhanced by the factor acting on the 
exponential leading contribution ([8]); this leads in particular to the rescaling of the anomalous dimension 
7med (7 — > 7med ~ \/^7o) cquivalcntly, to the enhancement of the in medium coupling constant. 
Since hard corrections are suppressed by the extra factor 1 / \/]V7, it is straightforward to check that r, G2 
and Q2 approach the asymptotic DLA limits tq = Nc/Cp, G2 = 4/3 and Q2 = ^ + Nc/3Cf [4] when 
Ns increases. The previously mentioned KNO-scaling experienced by G2 and Q2 proves no special 
sensibility to the model and should normally hold like in the vacuum. 

Finally, since these results are model-dependent, they may still be improved in the future, specially after 
the A'^s -dependence of the non-singular parts of the splitting functions (01) has been exactly computed. 

Perspective: Many experimental characterizations of the medium-modified intrajet structure in heavy- 
ion collisions at RHIC and at the LHC require a soft momentum cut-off pj^^, with Q > ' to remove the 
effects of the high multiplicity background. In [3], the soft background was subtracted by integrating the 
single inclusive differential distribution ^j^^ ("hump-backed plateau") over the range Q > pr ^ Pt^, 
with p™' > Aqcd- Accordingly, the equivalent computation should be performed for the second mul- 
tiplicity correlator by integrating the double differential inclusive distribution (two-particle correlation) 
dinpi'^rdinp2 T ^^^'^ "^^^^ lower bounds of integration p^"^ > Aqcd (^ = 1> 2). Imposing such 
a cut-off in our calculations will affect the normalization rather than the behavior and the shape of these 
observables as a function of Ns and the jet energy scale of the process Q [14]. 
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